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Abstract 

Using modern nucleon-nucleon interactions in the description of the ^4 = 3,4 nuclei, it is not 
possible to reproduce both the three- and four-nucleon binding energies simultaneously. This is 
one manifestation of the necessity of including a three-nucleon force in the nuclear Hamiltonian. In 
this paper we will perform a comparative study of some, widely used, three-nucleon force models. 
We will analyze their capability to describe the aforementioned binding energies as well as the n — d 
doublet scattering length. A correct description of these quantities can be considered a stringent 
requirement for a nuclear Hamiltonian containing two- and three-nucleon interaction terms. As 
we will show, this requirement is not fulfilled by several of the models available in the literature. 
To satisfy it, we propose modifications in the parametrization of the three-nucleon forces and we 
study their effects on few selected N — d low energy scattering observables. 

PACS numbers: 21.30.-x,21.45.Ff,27.10.+h 
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I. INTRODUCTION 



Realistic nucleon-nucleon (NN) potentials reproduce the experimental NN scattering data 
up to energies of 350 MeV with a x 2 P er datum close to 1. However, the use of these 
potentials in the description of the three- and four-nucleon bound and scattering states 
gives a x 2 P er datum much larger than 1 (see for example Ref. [l]). In order to improve that 
situation, different three-nucleon force (TNF) models have been introduced so far. Widely 
used in the literature are the Tucson-Melbourne (TM) and the Urbana IX (URIX) models 

n n 

[2|, |3J. These models are based on the exchange mechanism of two pions between three 
nucleons. The TM model has been revisited within a chiral symmetry approach [41], and it 
has been demonstrated that the contact term present in it should be dropped. This new 
TM potential, known as TM', has been subsequently readjusted [5]. The final operatorial 
structure coincides with that one given in the TNF of Brazil already derived many years 
ago [61. TNF models based on rrp and on meson exchange mechanisms have also been 

q n 

derived [7| and their effects have been studied in the triton binding energy |8(. More recently, 
TNFs have been derived using a chiral effective field theory at next-to-next-to-leading 
order. A local version of these interactions (hereafter referred as N2LOL) can be found 

n 

in Ref. [10]. At this particular order, the TNF has two unknown constants that have to 
be determined. More in general, all the models contain a certain number of parameters 
that fix the strength of the different terms that compose the interaction. It is a common 
practice to determine these parameters from the three- and four-nucleon binding energies. 
In the chiral effective field theory there is a consistent derivation of the two- and three- 
nucleon interactions and some of the low energy constants entering in the TNF are fixed 
already from the NN data. On the other hand, the parametrization of the TM' and URIX 
interactions have been determined in association with specific NN potentials. Therefore, 
their parametrizations could change when used with different NN potentials since different 
NN potentials predict different A = 3,4 binding energies. 

The n — d doublet scattering length 2 a n d can give valuable information. In principle 
this quantity is correlated, to some extent, to the A = 3 binding energy through the so- 
called Phillips line 11, Ijj]. However the presence of TNFs of the type studied here breaks 



this correlation. Therefore 2 a n d emerges as an independent observable that can be used to 
evaluate the capability of the interaction models to describe the low energy region. Due 
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to the lack of excited states in the A = 3 system, the zero energy state is the first one 
above the ground state. In the case of n — d scattering at zero energy, the J = \ + state 
is orthogonal to the triton ground state and, for this reason, the wave function presents a 
node in the relative distance between the incident nucleon and the deuteron. The position 
of the node is related to the scattering length and it is also sensitive to the relation between 
the overall attraction and repulsion of the interaction. Several of the realistic NN potentials 
underestimate the triton binding energy. Adding a TNF, which in general can include an 
attractive as well as a repulsive component, with a strength fixed for example to reproduce 
the triton binding energy, the balance between the overall attraction and repulsion of the 
interaction changes with respect to that one produced by the NN potential alone. And, as 
we will show, this leads to different predictions of 2 a n d and the a-particle binding energy 
i?( 4 He). An analysis of the parametrization of a chiral TNF, in order to describe the triton 
binding energy i?( 3 H), S( 4 He) and 2 a n d, has been performed in Ref. [9|. A similar analysis 
has not been done for the local models URIX, TM' and N2LOL since only the three- or 
four-body binding energy has been considered in the determination of their parametrization 
but not 2 (Ind- 



ia Ref. 



13] results for different combinations of NN interactions plus TNF models are 



given. We report the results for the quantities of interest in Table I. From the table, we 
can observe that the models are not able to describe simultaneously the A = 3,4 binding 
energies and 2 a n d- Triggered by this fact, in this paper we make a comparative study of the 
aforementioned TNF models. To this end we use the AV18 as the reference NN interac- 
tion and the three-nucleon interaction models will be added to it. Parametrizations of the 
URIX and TM' models already exist in conjunction with the AV18 potential. Conversely the 



N2LOL TNF has been constructed using the N3LO-Idaho potential from Ref. [15|. So, in a 
first step, we have adapted its parametrization in order to reproduce, in conjunction with the 
AV18 interaction, S( 3 H). Successively, we study the sensitivity of different parametrizations 
in the description of _B( 4 He) and 2 a n d- Selecting those parametrizations that predict these 
three quantities close to their experimental values, we study some polarization observables 
in p — d scattering at Ei a b = 3 MeV. As an interesting result, we have observed that the 
predictions of the different parametrizations fall in a narrow band that, in the case of the 
vector analyzing powers, has a different position for each model, indicating a sensitivity to 
the short range structure of the TNF. 
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All calculations have been done using the hyperspherical harmonics (HH) method as 



developed 
systems 



yy some of the authors to describe bound and scatte ring states in A = 3,4 



16 



19( | in configuration space or in momentum space 20|, |21J (for a recent review 



see Ref. 13|). The paper is organized as follows. In the next section we introduce the TNF 
models in configuration space defining their parametrizations. In Section III we make a 
sensitivity study of the parametrization for each model looking at _B( 3 H), B( 4 He) and 2 a n d- 
In Section IV we study p — d polarization observables at Ei a j, = 3 MeV for specific values of 
the parameters. The conclusions are given in the last section. 



II. THREE NUCLEON FORCE MODELS 



In Ref. [13j the description of bound states and zero-energy states for A = 3,4 has been 
reviewed in the context of the HH method. In Table [I] we report results for the triton and 
4 He binding energies as well as for the doublet n — d scattering length 2 a n d using the AV18 
and the N3LO-Idaho NN potentials and using the following combinations of two- and three- 
nucleon interactions: AV18+URIX, AV18+TM' and N3LO-Idaho+N2LOL. The results are 
compared to the experimental values also reported in the table. Worthy of notice is the 
recent very accurate datum for 2 a n d 22]. 

From the table we may observe that only the results obtained using an interaction model 
that includes a TNF are close to the corresponding experimental values. In the case of the 
AV18+TM', the strength of the TM' potential has been fixed to reproduce the 4 He binding 
energy and, as can be seen from the table, the triton binding energy is underpredicted. 
Conversely, the strength of the URIX potential has been fixed to reproduce the triton binding 
energy, giving too much binding for 4 He. The strength of the N2LOL potential has been 
fixed to reproduce simultaneously the triton and the 4 He binding energies. In the three cases 
the predictions for the doublet scattering length are not in agreement with the experimental 
value, in particular for the AV18+URIX model. 

Our intention is to study different parametrizations of the TNFs to obtain, as close as 
possible, a simultaneous description of the three quantities under observation. To this aim 
we give a brief description of the TM' (or Brazil), URIX and N2LOL models. Starting from 
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the following general TNF 

W= W (hJ,k) , (1) 

i<j<k 

a generic term can be put in the following form: 

W(l, 2, 3) = aW a (l, 2, 3) + bW b (l, 2, 3) + dW d (l, 2, 3) + c D W D (l, 2, 3) + c E W E (l, 2, 3) . (2) 

Each term corresponds to a different mechanism and has a different operatorial structure. 
The first three terms arise from the exchange of two pions between three nucleons. The 
a-term comes from irN S-wave scattering whereas the 6-term and d-term, which are the 
most important, come from nN P-wave scattering. The specific form of these three terms 
in configuration space is the following: 

W a {l, 2, 3) = Woin • r 2 )(cr 1 • r 31 )(<x 2 • r 23 )y(r 31 )y(r 23 ) 
W 6 (l,2,3) = W Q {n • t 2 )[(<7i • o-2)y(r 3 i)y(r 23 ) 

+ • r 3l )(rr 2 • r 23 )(r 3 i ■ r 23 )t(r 31 )t(r 23 ) 

+ (oi • r 31 )(cr 2 • r 31 )t(r 31 )y(r 23 ) 

+ (cri • r 23 )(cr 2 • r 23 )y(r 3 i)t(r 23 )] (3) 
W d (l,2,3) = W (t 3 ■ ri x r 2 )[(cr 3 • <x 2 x cTi)y(r 3l )y(r 23 ) 
+ • r 3 i)(o2 • r- 23 )(cr 3 • r 31 x r 23 )t(r 31 )t(r 23 ) 

+ (©"I ' r 3 i)(<T 2 • 1*31 X 0"3)i(^3l)y(^23) 
+ (0"2 • r 23 )(<T 3 • ^23 x 0"i)y (rai)* (r 23 )] , 

with W'o an overall strength. The b- and d-terms are present in the three models whereas 
the a-term is present in the TM' and N2LOL and not in UR1X. In the first two models, the 
radial functions y(r) and t(r) are obtained from the following function 

12* 1 [°° 2 jpfar) 
Mr)= m°wl dq " 7T^ FA{q) (4) 

where is the pion mass and 

{ (5) 

= ~.y'^ r ) ■ 

The cutoff function F A in the TM' or Brazil models is taken as F A = [(A 2 - m 2 )/(A 2 + g 2 )] 2 . 
In the N2LOL model it is taken as Fa = exp(— g 4 /A 4 ). The momentum cutoff A is a 



parameter of the model fixing the scale of the problem in momentum space. In the N2LOL, 
it has been fixed to A = 500 MeV, whereas in the TM' model the ratio A/m n has been 
varied to describe the triton or 4 He binding energy at fixed values of the constants a, b and 
d. In the literature the TM' potential has been used many times with typical values around 
A = 5 m n . 

In the URIX model the radial dependence of the b- and d-terms is given in terms of the 
functions 

Y(r) = e~ x /x^ Y 



T(r) = (1 + 3/x + 3/x 2 )Y(r) £ 



(6) 



T 



with x = m n r and the cutoff functions are defined as £y = £t = (1 — e~ cr ), with c = 2.1 
fm~ 2 . This regularization has been used in the AV18 potential as well. Since the URIX 
model has been constructed in conjunction with the AV18 potential, the use of the same 
regularization was a choice of consistency. The relation between the functions Y(r),T(r) 
and those of the previous models is: 

Y(r) = y(r) + T(r) 

(7) 

T (r) = -t(r) . 

With the definition given in Eq.(j4]), the asymptotic behavior of the functions fo(r), y{ r ) an d 
t(r) is: 

3 e~ x 

fo(r -> oo) -> 



ml x 



3e~ x ( 1 

y(r -)• oo) -> 1 + - ] (8) 



x 2 \ x 



3 e~ x / 3 3 



t(r -4ooM n H h 



nn2 rjQ V r^* 



To be noticed that with the normalization chosen for / , the functions Y and T defined from 

y and t and those ones defined in the URIX model coincide at large separation distances. 

Conversely, they have a different short range behavior. Using the URIX Y(r),T(r) functions, 

I I 

the a-term has been included in the construction of the Illinois TNF model [23] . 

The last two terms in Eq. (j2J) correspond to a two-nucleon (2N) contact term with a pion 
emitted or absorbed (D-term) and to a three-nucleon (3N) contact interaction (E-term). 
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Their local form, in configuration space, derived in Ref. [loj], is 

Wn{l,2,3) = W D (n ■ t 2 ){(<t 1 ■ <T 2 )[y(r 31 )Z (r 23 ) + Z (r 31 )y(r 23 )} 
+ (ai ■ r 3 i)(a 2 ■ r 31 )t(r 31 )Z (r 23 ) 

(9) 

+ (<ti ■ r 23 )(a 2 ■ r 23 )Z (r 31 )t(r 23 )} 
W E (1, 2, 3) = Wg{r x ■ r 2 )Z (r 31 )Z (r 23 ) . 

The constants W D and Wq fix the strength of these terms. In the case of the URIX model 
the D-term is absent whereas the E'-term is present without the isospin operatorial structure 
and it has been included as purely phenomeno logical, without justifying its form from a 
particular exchange mechanism. Its radial dependence has been taken as Z (r) = T 2 (r). In 
the N2LOL model, the function Zq(t) is defined as 

127T 1 /*°° 

Mr) = — / dqq 2 3o (qr)F A (q) (10) 

m 7T ATC JO 

with the same cutoff function used before, F\(q) = exp(— g 4 /A 4 ). In the TM' model the D- 
and E-terms are absent. 

Each model is now identified from the values assigned to the different constants. Following 
Refs. , in the case of the TM' model, the values of the constants are a = —0.87 m w , 

b = —2.58 m~ 3 , and d = —0.753 m~ 3 ; the strength Wo = {gm^/^m^) 2 m\ and the cutoff 
has been fixed to A = 4.756 m n in order to describe correctly, associated to AV18, I?( 4 He). 
In Tabled] the calculations have been done using these values with g 2 = 197.7, = 139.6 
MeV, m^/m^ = 6.726 (m^ is the nucleon mass) as given in the original derivation of the 
TM potential. 

In the URIX model the b- and ci-terms are present, however with a fixed ratio based on 
the Fujita-Miyazawa diagram. The strength of these terms are: bWo = 4 A™ and d = 6/4, 
with A™ = —0.0293 MeV. The model includes a purely central repulsive term introduced 
to compensate the attraction of the previous term, which by itself would produce a large 
overbinding in infinite nuclear matter. It is defined as 

WE mx {l, 2, 3) = A R T 2 {r 31 )T 2 {r 23 ) (11) 

with A R = 0.0048 MeV. 

In the N2LOL potential the constants of the a-, b-, d-, D- and E-terms are defined in the 
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following way: 



2 2 
9A m ir 
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w 

< = ^fe) (?)^ (12) 



12tt 2 \F 7 

1 / 77 

with a = cim 2 , 6 = c 3 /2, d = c 4 /4, and cj = -0.00081 MeV" 1 , c 3 = -0.0032 MeV' 1 , 
c 4 = —0.0054 MeV -1 taken from Ref. |l5l| . The other two constants, cd = 1.0 and eg = 
—0.029, have been determined in Ref. [10| from a fit to 5( 3 H) and -B( 4 He) using the N3LO- 
Idaho+N2LOL potential model. The numerical values of the constant entering in Wo, W D 
and W E are m n = 138 MeV, = 92.4 MeV, qa = 1.29, and the chiral symmetry breaking 
scale A x = 700 MeV. 

In order to analyze the different short range structure of the TNF models, in Fig. [T] 
we compare the dimensionless functions Zo(r), y(r) and T(r) for the three models under 
consideration. In the TM' model using the definition of Eq. ffTUl) and using the corresponding 
cutoff function we can define: 



TM 12tt 1 r° f ^-ml \ 2 3/m^/A 2 \ 2 _ Ar 

z ° (r) = ^2^ y dqq 3 ° {qr) (ah7 j = 2 ixH< " J • (l3) 

This function is showed in the first panel of Fig. [T]as a dashed line. From the figure we can 
see that, in the case of the URIX model, the functions Z (r) and y(r) go to zero as r — > 0. 
This is not the case for the other two models and is a consequence of the regularization 
choice of the Y and T functions adopted in the URIX. 



III. PARAMETRIZATION STUDY OF THE THREE NUCLEON FORCES 

In this section we study possible variations to the parametrization of the TNF models in 
order to describe the A = 3, 4 binding energies and 2 a n( i. 



A. Tucson-Melbourne Force 



We first study the TM' potential and we would like to see whether, using the AV18+TM' 
interaction, it is possible to reproduce simultaneously the triton binding energy and the 
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doublet n — d scattering length for some values of the parameters. The a-term gives a very 
small contribution to these quantities, therefore, in the following analysis we maintain it 
fixed at the value a = —0.87 m^ 1 . The analysis is shown in Fig. [2J In the left panel, the 
doublet n — d scattering length, 2 a nc i, is given as a function of the parameter b (in units of 
its original value bo = —2.58 m~ 3 ) for different values of the cutoff A (in units of m n ). The 
box in the figure includes those values of 2 a n d compatible with the experimental results. On 
each point of the curves, the value of the constant d has been varied to reproduce the triton 
binding energy. Its corresponding values (in units of its original value do = —0.753 m~ 3 ) are 
given in the right panel as a function of b. Therefore, each point of the curves in both panels 
corresponds to a set of parameters that, in connection with the AV18 potential, reproduces 
the triton binding energy. The variations of the parameters given in Fig. [2] do not exhaust all 
the possibilities. The analysis has been done maintaining the attractive character of the b- 
and d-terms and, therefore, the lines in the left panels of the figure stop when one of the two 
parameters, b or d changes sign. We can observe that, with the AV18+TM' potential, there 
is a very small region in the parameters' phase space available for a simultaneous description 
of the triton binding energy and the doublet scattering length. This small region corresponds 
to a value of b around four times bigger than the original value &o an d d results to be almost 
zero. Moreover, the value of the cutoff A around 3.8m„- is smaller than the values usually 
used with the TM' potential (A « 5m n ). 

To be noticed that, for negative values of the parameters a, b and d, the TM' potential 
is attractive and it does not include explicitly a repulsive term. Added to a specific NN 
potential that underestimates the three-nucleon binding energy, it supplies the extra binding 
by fixing appropriately its strength. As mentioned in Sec. I, the scattering length is sensitive 
to the balance between the attractive part and the repulsive part of the complete interaction. 
Therefore, in the case of the TM' potential, it seems that introducing only attractive terms, 
fixed to reproduce the triton binding energy, it is difficult to reproduce correctly this balance. 

As discussed before, the TM' potential is a modification of the original TM potential 



compatible with chiral symmetry. At next-to-next- 
field theory the D- and E-terms appear (see Ref. 



;o-leading order in the chiral effective 
9J and references therein) as given in 



Eq.([2"l). Here we introduce the following additional term to the TM' potential based on a 



9 



contact term of three nucleons 

W™{1, 2, 3) = W E Z™(r 31 )Z™(r 23 ) . (14) 

This term corresponds to the i?-term in Eq. (j2J), except that, for the sake of simplicity, we 
have omitted the (t\ ■ Tq) operator. Its strength Wq is defined in Eq. ffl2|) and the function 
Zq M , defined in Eq. ( fT3|) . is a positive function, therefore, for positive values of ce, the new 
term is repulsive. We include it in the following analysis of the TM' potential. The results 
are shown in Fig. [3] for three values of A/m n = 4, 4.8, 5.6. In the left panels the doublet n — d 
scattering length is given as a function of the parameter b (in units of b ) for different values 
of the strength ce of the i?-term. The box in the panels includes those values compatible 
with the experimental results. At each point of the curves, the value of the constant d has 
been varied to reproduce the triton binding energy. For selected values of the parameters 
inside the box, the predictions for the 4 He binding energy, i?( 4 He), are shown in the right 
panels. 

Comparing the left panels in Figs. [2] and [3], the effect of the new term is clear. In Fig. [2] we 
have observed that using A > 4 m n , 2 a n d cannot be well reproduced. Conversely, in Fig. [3l 
the inclusion of the new term allows for a description of 2 a n d using different values of the 
cutoff. The values of the parameter b are closer to its original value as A increases. Opposite 
to this, the predictions of i?( 4 He) improves as A decreases. For example, considering the 
case A = 4 m n , 2 a nd , £>( 3 H) and I?( 4 He) are well reproduced with b = 3.26 , d = 6.2d and 
Ce = 1- With A = 4.8 the set of parameters that gives the best description of the three 
quantities is b = 1.5&o, d = 4.5<io and ce = 1-6. And with A = 5.6 they are b = 0.8b , 
d = 3do and ce = 2. Their different contributions to the triton binding energy are given 
in Table [III where we report the mean values of the kinetic energy and the NN potential 
energy as well as the mean values of the attractive part of the TNF, Va{3N), corresponding 
to the sum of the a, b and d-terms, and the repulsive part, Vr(3N), corresponding to the ce 
term. The last two columns show i?( 4 He) and 2 a n d- For the sake of comparison, in the first 
row, the original values of the parameters have been considered (b = bo,d = d and eg = 0) 
with the value of the cutoff fixed to reproduce the triton binding energy (A = 4.8 m^). As 
we can observe, in this case _B( 4 He) is overestimated and 2 a n d is underestimated. When the 
E-term is considered, the description of _B( 4 He) improves and it seems that a low value of 
A is preferable. A further analysis of these parametrizations is given Sec. IV studying some 
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polarization observables at low energy. 
B. Urbana IX Force 

In the following we analyze the URIX potential which has two parameters, called A™ 
and Ar. In this model the strength of the ci-term is related to the strength of the 6-term as 

n 

d = b/A. The original values of the parameters have been fixed in Ref. [3j in conjunction with 
the AVI 8 NN potential and, from Table HJ we observe that the model correctly describes 
the triton binding energy. However, it overestimates _B( 4 He) and underestimates 2 a n d- In 
order to further analyze the origin of this behavior, we have varied the constants A™ , Ar 
and the relative strength = d/b of the b- and <i-terms. The regularization parameter 

has been held fixed at its original value, c = 2.1 fm -2 . For a given value of A™ , the values 
of Ar and has been chosen to reproduce S( 3 H) and 2 a n d- The results are shown in 

Fig. HI In panel (a), is given as a function of A™ with Ar varying from 0.0176 MeV 

at A™ = -0.02 MeV to 0.0210 MeV at A%™ = -0.050 MeV. These values of A R are 
more than three times bigger than the original value of 0.0048 MeV. In panel (b) and (c) 
the results for 2 a n d and _B( 4 He) are given respectively. The latter has not been included in 
the determination of the parameters, since and Ar have been determined from the 

triton binding energy and 2 a n d, and is therefore a pure prediction. We observe a slightly 
overestimation of -B( 4 He), in particular for values of l^fJH > 0.04 MeV, corresponding to 
values of D™ < 0.7. 

With modifications of the parameters in the URIX force, we were able to describe rea- 
sonably well .B( 3 n), 2 a n d and -B( 4 He). However, this has been achieved with a substantial 
increase of the repulsive term. In order to gain insight on the consequence of the new 
parametrizations in the quantities of interest, in Table [TTT1 we report the mean values of the 
kinetic energy and the NN potential energy as well as the mean values of the attractive part 
of the TNF, Va{3N), corresponding to the sum of the b and d-terms, and the repulsive part, 
Vr(3N), corresponding to the Ar term, for selected values of the parameters (indicated as 
points in Fig. H|). The last two columns show I?( 4 He) and 2 a n d- For the sake of comparison, 
in the first row, the values obtained using the original AV18+URIX model are reported. 
From the table we observe that some of the values considered for an d Ar are quite far 
from the original ones. At the original value of A™, —0.0293 MeV, the relative strength 
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now results to be D^y = 1 and Ar = 0.0181 MeV. As diminishes, Ar tends to in- 

crease further with the consequence that the mean value Vr(3N) is more than three times 
larger than the value obtained using the original parameters (given in the first row). This is 
compensated by a lower mean value of the kinetic energy. A further analysis of the effects 
of the parametrizations given in Table IIHI is performed in Sec. IV studying selected p — d 
polarization observables. 



C. N2LOL Force 



The parameters ci, C3 and C4 of the N2LOL model have been taken from the the chiral 
N3LO NN force of Ref. 15], whereas the cd and ce parameters have been determined in 



Ref. 



10], in conjunction with that NN force, by fitting B( 3 H) and _B( 4 He). Here we are 



going to use the N2LOL force in conjunction with the AV18 NN interaction, so we have 
to modify its parametrization since the amount of attraction to be gained is now different 
(see Tabled]). In the following we will call c°, c° and c° A the values of these constants, given 



in Sec. II, determined in Ref. 



15] . Among different possibilities, in Fig. [5] we show a new 



parametrization of the N2LOL interaction obtained by multiplying c° and c° by a factor cq 
and maintaining c\ = c\. Then the parameters Cd and ce have been determined from a fit 
to 5( 3 H) and 2 a nd . They are shown on panel ( function of cq. Therefore, at a fixed 

value of Co, with the set of parameters c\ = c°, C3 = cqc\, C4 = c$c° 4 and the corresponding 
values of cp and ce extracted from the figure, the AV18+N2LOL interaction reproduces 
the B( 3 H) and 2 a n d- In panel (b) we show the stability obtained in the description of the 
doublet scattering length corresponding to the constant value chosen for the determination 
of the parameters, 2 a n d = 0.644 fm. With the new set of parameters it is now possible to 
calculate _B( 4 He). This is shown in panel (c) and it is interesting to note that in all cases the 
value .B( 4 He)=28.60 MeV has been obtained. Modifying also the parameter c\ as c% = cqc\ 
slightly different values of cd and ce are obtained. Using these values to calculate -B( 4 He), 
again we obtain a constant value that is now = 28.55 MeV. Similar analyses using slightly 
different values of the cutoff A around 500 MeV do not change these results. 

In order to correctly describe _B( 4 He), after fixing i?( 3 H) and 2 a n d, we now analyze a 
modification to the relative strength of the b- and d-terms which, in the previous analysis, 
was maintained at its original value of C4/C3 = 1.6875. To this end we perform a similar 
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study as has been done previously for the other TNF models. Fixing the constant c\ to its 
original value c°, C3, C4, cd and ce have been varied. The analysis is shown in Fig. [6] at the 
following four values of ce = 0,0.1, —0.03, —0.5 and for the indicated values of cp chosen 
to reproduce 2 a n d (left panels). The predictions for S( 4 He) are given in the right panels 
for those values of the parameters that give a value of 2 a n d inside the box. At each point 
of the curves in the left panels and at the points in the right panels, C4 has been chosen 
to reproduce the triton binding energy. Due to the (r^ ■ Tj) operator in the i?-term of the 
N2LOL potential, positive values of ce makes this term attractive. Conversely, negative 
values of ce makes this term repulsive. We have considered only one positive case, ce = 0.1. 
Increasing further ce we found it difficult to describe correctly _B( 4 He). For negatives values 
of Ce we have considered two cases, ce = —0.03, which corresponds to the value given in 
Ref. [loj], and ce = —0.5. From the figure we observe an almost linear behavior of 2 a n d- 
There is a slight curvature for negative values of cd in the upper three panels. The analysis 
of _B( 4 He) selects the values of cb- We have found that the experimental value is well 
reproduced for the pairs (cd = — 0.5, ce = 0.1), (cd = — 1,ce = 0), (cd = — 1,ce = —0.03), 
and (cd = — 2,cg = —0.5). 

In Table IIVI we report the mean values of the kinetic energy, the two-nucleon potential 
energy as well as the mean values of the attractive part of the TNF, Va[3N), and its 
repulsive part Vr(3N) for the selected values of the parameters that correspond to the best 
description of the three quantities under study. In the last two columns of the table, _B( 4 He) 
and 2 a n d are given. The contributions to Va{3N) come from the a-, b-, d and D-terms, which 
are always attractive in the cases considered, and from the E-teim. in the first case. This 
term contributes to the repulsive part Vr(3N) in the last two cases. From the table we may 
observe that C3 and C4 results to be larger and smaller than their original values, respectively. 
This is a consequence of the simultaneous description of S( 3 H) and 2 a n d- Furthermore, in 
the first three cases, the ratio C4/C3 ~ 0.46, is much smaller than the original ratio. 

IV. ANALYSIS OF THE POLARIZATION OBSERVABLES 

In the previous section we have studied different parametrizations of the TM', URIX 
and N2LOL TNF models in conjunction with the AV18 NN potential. The analysis has 
been done varying the parameters in order first to reproduce B( 3 H) and then looking at 
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their dependence on 2 a n( i and £>( 4 He). To improve the description of these quantities, some 
substantial modifications were necessary for the first two models. In the case of the TM' 
interaction we found opportune the inclusion of a repulsive term. In the analysis of the URIX 
interaction, the strength of the repulsive term resulted to be more than three times bigger 
than the original value and the relative strength of the b- and <i-terms, originally fixed to 
1/4, has been also increased. In the case of the N2LOL interaction, some adjustment of the 
parameters was necessary, mainly due to the fact that the AV18 interaction is less attractive 
than the N3LO interaction, from which the N2LOL model has been originally parametrized. 
In this section we analyze the effects of the new parametrizations in observables that are 
not correlated to the binding energies or to 2 a n ^. Some polarization observables in p — d 
scattering have this characteristic, in particular the vector and tensor analyzing powers. In 
Figs. [71 [U |9] we show the differential cross section da/dfl, the vector polarization observables 
A y and iT\\ and the tensor polarization observables T 2 o, T 2 \ and T 22 at Ei a f, = 3 MeV for the 
different potential models compared to the results obtained using the original AV18+URIX 
interaction. In the figures, the cyan band collects the results obtained with the parameters 
given in the three last rows of Table [Til from the second to the sith row of Table IIII1 and the 
three last rows of Table ITVl for each model respectively, whereas the solid line is the prediction 
of the original AV18+URIX model. As we can see, for each TNF model, the observables 
calculated using the different parametrizations, fixed from a simultaneous description of 
5( 3 H), 2 a n d and _B( 4 He), fall in bands which, in the case of the vector analyzing powers, 
have a different position for the three models. Since the models essentially differ in the 
definitions of the functions y(r), T(r) and Zo(r), this difference can be associated to the 
different short-range behavior of the TNF models. In Fig. [7]we observe that the AV18+TM' 
model, using the new parametrizations, does not give any improvement in the observables 
compared to the AV18+URIX predictions. Moreover iT n and T 2 i are worse described. It 
should be observed that the AV18+TM' model, with the original parametrization, and the 
AV18+URIX give similar results for the observables (a small difference can be observed in 
the maximum of A y being slightly higher for the former). Therefore the previous conclusions 
do not change if compared to the original AV18+TM' model. In Fig. [8] we observe that the 
new parametrizations of the AV18+URIX produce a much worse description of A y , iT n and 
T 2 i. Since the vector analyzing powers are mainly described by the P-wave phase-shift and 
mixing parameters, we can conclude that they result to be poorly reproduced with the new 
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parametrizations. Conversely to what happened analyzing the previous models, in Fig. [9], 
we observe that the N2LOL interaction produces an improvement in the description of A y 
and iT\\. The well known discrepancy in these observables is now reduced and, in particular 
for A y , the improvement is noticeable. In the case of the tensor analyzing powers, a slightly 
worse description of T 2 i between the two maxima is now observed. In general all TNFs of the 
type analyzed here have this effect in T21 indicating that a different mechanism, not present 
in the models, should be considered to improve the description of the minimum around 75°. 

Finally we would like to comment on the fact that the vector analyzing powers, A y and 
iT\\, calculated using different TNF models fall inside a band with a different position for 
each model. In Fig. (TUlthe three bands, extracted from Figs. IT1I81I91 are shown explicitly and 
compared to original AV18+URIX model (solid line). We can clearly observe the different 
position of the bands with the best description obtained with the new parametrizations 
of the AV18+N2LO model and the worst description with those one of the AV18+URIX 
model. Since all the models inside the bands describe reasonably well £>( 3 H), 2 a n d and 
_B( 4 He), we can conclude that the difference is a direct consequence of their different short 
range structure. A natural question is whether, with opportune modifications of their radial 
dependence, i.e. modifying the functions y(r), T(r) and Z (r), it will be possible to improve 
further the description of these observables at Euj, = 3 MeV and, eventually, obtain a \ 2 
per datum close to one. A preliminary study in this direction has shown that a further 
improvement in the A y and %T\\ maxima is associated to a worse description of the T 2 \ 
minimum. The particular structure of these observables is related to a bigger splitting 
in the 4 Pj phase-shifts than the normal splitting produced by the two-nucleon forces, as 
26| . In particular the 4 Pi/2 phase-shift has to be smaller and the mixing 



discussed in Ref. 



parameter €3/2- has to be bigger. It is a general feature of the TNFs studied here that 
they tend to increase both, A P\/2 and e 3 / 2 -. To be more precise, in Table [V] we show 
the A Pj phase-shifts and e 3 / 2 - for the AV18 and AV18+URIX potential models, and for one 
selected set of the parameters of Tables HT1 IHIl HVl corresponding to the new parametrizations 
of the AV18+TM', AV18+URIX and AV18+N2LOL models (indicated in the table with 
an asterisk). In particular, parametrizations of the second row of Table [Til fourth row of 
Table [TTTJ and third row of Table HVl have been used, respectively. In the last row of the table, 



the results from phase-shift analysis (PSA) of Ref. 26[ are given. From the table we observe 
that the 4 Pi/2 phase-shift increases when the TNF models are added to the AV18 potential. 
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By itself this change will produce a much worse description of A y and %T\\. However this 
is well compensated with the corresponding increase in 4 Ps/2 and e 3 / 2 -. This is not the 
case with the minimum in T 2 i, for which a better description would be obtained lowering 
the AV18 value of 4 Pi/2, as discussed in Ref. 26]. The other parametrizations given in 
Tables [TTJ 1111} IIVI produce similar changes in the 4 Pj parameters. From this observation we 
can conclude that the spin-isospin structure of the TNFs considered here is not sufficient to 
describe simultaneously _B( 3 H), 2 a n d, -B( 4 He) and the vector and tensor analyzing powers at 
low energies. 



V. CONCLUSIONS 



Stimulated by the fact that some of the widely used TNF models do not reproduce 
simultaneously the triton and the 4 He binding energies and the n — d doublet scattering 
length, we have analyzed possible modifications to their parametrizations. To this end we 
have selected the AV18 as the reference two-nucleon force and, associated with it, we have 
varied the original parameters of the TM' and URIX models so as to improve the description 
of the three quantities mentioned. Furthermore, using the recent local form of a chiral TNF 
(we have called this model N2LOL), we have studied its parametrization associated to the 
AV18 interaction too. The analysis has proceeded in the following way. The three models 
under observation, TM', AV18 and N2LOL, have been written in configuration space as a 
sum of five terms, the a-, b-, d-, C- and P-terms. The first three, corresponding to a two- 
pion exchange process, are attractive. The last two, corresponding to contact terms, can be 
either attractive or repulsive. Not all the models include the five terms. In the TM' model 
only the a-, b- and d-terms are present and therefore, this model does not include explicitly 
a repulsive term. The URIX model includes the b-, d- and P-terms. This last term has 
been parametrized as repulsive in order to compensate the large overbinding produced by 
the first two terms in infinite nuclear matter. The N2LOL model includes the five terms. 

The study has been started analyzing the AV18+TM' model. Maintaining fixed the 
strength of the a-term at its original value, we have varied the strengths of the b- and d- 
terms for several values of the cutoff parameter A. We have explored negative values of the 
strength parameters b and d in order to keep the attractive character of these terms. We 
have found it difficult to reproduce 2 a n d for reasonable values of the strength parameters. 
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This fact has motivated the subsequent step of introducing a repulsive term in the model. 
As a simple choice, we have introduced a purely central E-term and a corresponding Zo(r) 
function, obtained using the monopole cutoff of the model. Including this term we were able 
to describe simultaneously B( 3 H) and 2 a n d for several values of the cutoff. A further selection 
among these values has been done from the calculation of -B( 4 He). We have observed that 
with A < 4.8 m n it was possible to describe the three quantities reasonably well. 

In the original AV18+URIX model the relative strength between the b- and d-terms 
was fixed. In the present analysis we have relaxed this condition increasing the number of 
parameter of the model from two to three, the strengths of the b-, d- and i?-terms. Varying 
them, we have found it possible to describe the three quantities of interest for values of 
the parameters very different from their original ones. In particular, the strength of the 
repulsive term resulted more than three times larger than the original value. In the case of 
the AV18+N2LOL model, maintaining the strength of the a-term fixed to its original value, 
we have varied the parameters C3, C4, cd and ce in combinations that reproduce 5( 3 H). 
Then we have studied the dependence on 2 a n d and B( A He) of the different parametrizations. 
For fixed values of eg we have calculated 2 a n d for different values of C3 and cd- We have 
found that C3 > 1.40° in order to describe simultaneously B( S H) and 2 a n d- The values of 
cd has been selected from the analysis of _B( 4 He). Values of -B( 4 He) compatible with the 
experimental value have been found in the four cases of eg explored. 

After making this sensitivity study we have selected, for each model, some combinations 
of the parameters that give the better description of B( 3 H), 2 a n d and -B( 4 He) and we have 
calculated the differential cross section and the vector and tensor analyzing powers at = 
3 MeV. At this energy there are well established discrepancies between the predictions 
of the theoretical models and the experimental results. For example all potential models 
underestimate A y (the so-called A y puzzle) and iTn and overestimate the central minimum 
in T21. Some TNF models have been constructed ad hoc to improve the description of 
these observables at low energy 27]] . However the models studied here, derived from the 
exchange of two pions and contact terms, are not able to solve these discrepancies. What we 
have observed in the present study is that after fixing the parameters of each model from the 
description of B( 3 H), 2 a n d and -B( 4 He), the description of the vector polarization observables 
lies in a narrow band, positioned differently for each model. The best description is given by 
the AV18+N2LOL model which, with respect to the original AV18+URIX model, reduces 



17 



appreciably the discrepancy in A y and %T\\. However it gives a slightly worse description of 
the central minimum of Ti\. The other two models do not improve the description of the 
observables, compared always to the original AV18+URIX. The modified TM' model gives 
similar results though iT\\ is slightly worse whereas the results with the modified URIX 
model are definitely worse than the original model. The fact that for each model the A y 
and iTn predictions lie in a narrow band, indicates a connection between the short-range 
structure of the TNF and the polarization observables at low energies. From the analysis we 
can conclude that the smoother form of the y(r), T(r) and Z (r) functions of the N2LOL 
potential are preferable. To be noticed that the TM' and URIX models do not include a D- 
term. An extended analysis of these two models including it, will allow for a more stringent 
conclusion about the short-range structure of the TNF models. Preliminary studies in this 
direction are underway. 

Finally, at the end of Sec. Ill, we have analyzed the 4 Pj phase-shift parameters. The 
overall attractive character of the TNF makes larger the A P\/2 and 4 ^5/2 parameters, com- 
pared to the ones obtained using a NN force alone, and has little effect on 4 Ps/2- The mixing 
parameter €3/2- is larger too. Depending on the relative increase of these parameters, the 
description of A y and iT u can improve, as in the case of the N2LOL model, but not the 
description of the central minimum of T21. The spin-isospin structure of the TNF models 
studied here cannot lower the 4 Pi/2 phase-shift which seems to be necessary in order to 
improve the description of T21 in the minimum. A different mechanism has to be included 
in the structure of the TNF, as for example it has been proposed in Ref. 2^]. Studies along 
this line are at present underway. 
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TABLE I: The triton and 4 He binding energies B (MeV), and doublet scattering length 2 a n d (fm) 
calculated using the AV18 and the N3LO-Idaho two-nucleon potentials, and the AV18+URIX, 
AV18+TM' and N3LO-Idaho+N2LOL two- and three-nucleon interactions. The experimental 
values are given in the last row. 



Potential S( 3 H) £( 4 He) 2 a nd 



AV18 


7.624 


24.22 


1.258 


N3LO-Idaho 


7.854 


25.38 


1.100 


AV18+TM' 


8.440 


28.31 


0.623 


AV18+URIX 


8.479 


28.48 


0.578 


N3LO-Idaho+N2LOL 


8.474 


28.37 


0.675 


Exp. 


8.482 


28.30 


0.645±0.003±0.007 



TABLE II: Mean values of the triton kinetic energy and the two-nucleon potential energy V(2N), 
and the attractive, Va(3N), and repulsive, Vr(3N), contributions of the TNF to the triton binding 
energy using the AV18+TM' potential for the specified values of the parameters and with a = 
—0.87 m~ x . In the last two columns B{ He) and 2 a n d are given respectively. The experimental 
values are given in the last row. 



b 


d c E 


A 


T 


V(2N) V A {3N) V R (3N) £( 4 He) 


O-nd 








[MeV] 


[MeV] [MeV] 


[MeV] 


[MeV] 


[fm] 


-2.580 


-0.753 0.0 


4.8 


50.708 


-58.144 -1.039 


0.0 


28.52 


0.596 


-8.256 


-4.690 1.0 


4.0 


50.317 


-57.366 -2.206 


0.781 


28.30 


0.644 


-3.870 


-3.375 1.6 


4.8 


50.699 


-57.641 -2.748 


1.215 


28.38 


0.644 


-2.064 


-2.279 2.0 


5.6 


50.998 


-57.940 -2.814 


1.291 


28.44 


0.640 



Exp. 28.30 0.645±0.003±0.007 
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TABLE III: Mean values of the triton kinetic energy, the twonucleon potential energy V(2N), and 
the attractive, Va(3N), and repulsive, Vft(3iV), contributions of the TNF to the triton binding 
energy using the AV18+URIX potential for the specified values of the parameters. In the last two 
columns i?( 4 He) and 2 a n d are given respectively. The experimental values are given in the last row. 



A PW 


7^1 PW 


Ar T 


V(2N) V A (3N) V R {3N) 5( 4 He) 


2 

O-nd 


[MeV] 




[MeV] [MeV] 


[Me v j 


|lviC V J 


i ivie v j 


[ivie v j 


[fm] 


-0.0293 


0.25 


0.0048 51.259 


-58.606 


-1.126 


1.000 


28.48 


0.578 


-0.0200 


1.625 


0.0176 47.472 


-57.976 


-0.923 


2.950 


28.33 


0.644 


-0.0250 


1.25 


0.0182 47.628 


-57.967 


-1.162 


3.024 


28.34 


0.644 


-0.0293 


1.00 


0.0181 47.876 


-58.000 


-1.369 


3.015 


28.33 


0.643 


-0.0350 


0.8125 


0.0191 47.998 


-57.975 


-1.649 


3.147 


28.33 


0.645 


-0.0400 


0.6875 


0.0198 48.133 


-57.964 


-1.897 


3.249 


28.38 


0.645 


-0.0450 


0.5625 


0.0198 48.414 


-57.995 


-2.148 


3.248 


28.38 


0.643 


-0.0500 


0.50 


0.0210 48.471 


-57.952 


-2.401 


3.401 


28.44 


0.645 


Exp. 












28.30 


0.645±0.003±0.007 



TABLE IV: Mean values of the triton kinetic energy, the two-nucleon potential energy V(2N), and 
the attractive, Va(3JV), and repulsive, Vr(3N), TNF contributions to the triton potential energy 
using the AV18+N2LOL potential for the specified values of the parameters and with c\ = —0.00081 
MeV -1 . In the last two columns i?( 4 He) and 2 a n d are given. The experimental values are given in 
the last row. 



C3 


C4 


CD 


CE T 


V{2N) V A (3N) Vr(3N) B( 4 He) 


a n d 


[c§] 


[4] 




[MeV] 


[MeV] [MeV] 


[MeV] 


[MeV] 


[fm] 


1.4 


0.3636 


-0.5 


0.1 49.834 


-57.278 -1.029 


0.0 


28.31 


0.641 


1.4 


0.3786 


-1 


0.0 49.950 


-57.401 -1.022 


0.0 


28.30 


0.636 


1.5 


0.3735 


-1 


-0.03 49.839 


-57.274 -1.076 


0.036 


28.29 


0.644 


1.7 


0.9000 


-2 


-0.50 50.166 


-57.181 -2.119 


0.657 


28.32 


0.645 


Exp. 












28.30 


0.645±0.003±0.007 
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TABLE V: The 4 Pj phase-shifts and the £3/2- mixing parameter at Ei a b = 3 MeV for the potential 
models indicated. For the sake of comparison, the results of the PSA from Ref. {^(J are given in 
the last row. 





4 D 
Pl/2 


4 D 
"3/2 


4 D 
P5/2 


e 3/2- 


AV18 


22.03 


24.24 


24.08 


-2.247 


AV18+URIX 


22.31 


24.30 


24.27 


-2.314 


AV18+TM'* 


22.79 


24.45 


24.53 


-2.453 


AV18+URIX* 


22.75 


24.41 


24.35 


-2.375 


AV18+N2LOL* 


22.55 


24.25 


24.48 


-2.394 


PSA 


21.77 ±0.01 


24.30 ±0.01 


24.26 ±0.01 


-2.46 ±0.01 




FIG. 1: The Zq(t), y(r) and T{r) functions as functions of the interparticle distance r for the 
URIX (solid line), TM' (dashed line) and N2LOL (dotted line) models. 
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FIG. 2: The doublet scattering length 2 a n d as a function of the parameter b (in units of the 
original parameter bo = —2.58 m~ 3 ) of the TM' potential for different values of the cutoff and the 
corresponding values of the parameter d (in units of the original parameter do = —0.753 m~ 3 ), as 
a function of the parameter b, used to reproduce the triton binding energy. 
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FIG. 3: The doublet scattering length a n d as a function of the parameter b (in units of bo = 
—2.58 m~ 3 ) of the TM' potential including the W™ -term, for different values of the strength eg 
and for three selected values of A. The corresponding values of i?( 4 He), for specific values of a nc i 
inside the box, are shown in the right panels (circles). 
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FIG. 4: (a) The relative strength Dgj? , (b) 2 a nd and (c) 5( 4 He) as functions of A%™ , for the 
AV18+URIX model. 




FIG. 5: (a) The cd and ce parameters, (b) 2 a n d and (c) B( 4 He), as functions of cq, for the 
AV18+N2LOL model. 
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FIG. 6: The doublet scattering length 2 a n d as a function of the parameter C3 (in units of C3 = 
—0.0032 MeV -1 ) of the N2LOL potential, for different values of the strength cd, at four selected 
values of ce- The corresponding values of _B( 4 He), for specific values of a n d inside the box, are 
shown in the right panels (circles). 
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FIG. 7: (Color on line) Differential cross section and vector and tensor polarization observables 
at Eiat, = 3 MeV using the AV18+TM' model with the parameters given in the last three rows of 
Table ITT1 (cyan band). The predictions of the AV18+URIX model (solid line) and the experimental 
points from Ref. [25j are also shown. 
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FIG. 8: (Color on line) Differential cross section and vector and tensor polarization observables at 
Eiab = 3 MeV using the AV18+URIX model with the parameters given in Table Hill (cyan band). 



The predictions of the original AV18+URIX mode 
as a solid line. The experimental points from Ref. 



, given in the first row of the table, are shown 



251 ] are also shown. 



28 




FIG. 9: (Color on line) Differential cross section and vector and tensor polarization observables at 
Elab = 3 MeV using the AV18+N2LOL model with the parameters given in Table [IV] (cyan band). 



The predictions of the AV18+URIX model (solid line) and the experimental points from Ref. 
are also shown. 
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K 1(1. 10: (Color on line) The vector analyzing powers A y and iTu at Ei a b = 3 MeV using the 
AV18+TM' (cyan band), AV18+URIX (violet band) and AV18+N2LOL (red band) models as in 
Figs. I7I8I9| T he predictions of the AV18+URIX model (solid line) and the experimental points 



from Ref. 



25[ are also shown. 
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